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Abstract 

Using first-principles evolutionary simulations, we have systematically investigated phase sta¬ 
bility in the Hf-0 system at pressure up to 120 GPa. New compounds Hf 502 , Hf 302 , HfO and 
HfOa are discovered to be thermodynamically stable at certain pressure ranges and a new sta¬ 
ble high-pressure phase is found for Hf20 with space group Pnnm and anti-CaCl 2 -type structure. 
Both P62m-HfO and P 4 m 2 -Hf 203 show semimetallic character. Pnnm-HfOs shows interesting 
structure, simultaneously containing oxide and peroxide [0-0]^“ anions. Remarkably, it is 
P62m-HfO rather than 0II-Hf02 that exhibits the highest mechanical characteristics among Hf-0 
compounds. Pnnm-Hf 20 , /mm 2 -Hf 502 , P31m-Hf20 and P 4 m 2 -Hf 203 phases also show superior 
mechanical properties, these phases can be quenched to ambient pressure and their properties can 
be exploited. 
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I. INTRODUCTION 


Hafnium oxide Hf02 has a wide range of technological applications. In electronics indus¬ 
try, hafnium oxide-based material is currently used as an excellent high-/j gate dielectric^ 
and oxygen-deficient hafnium oxide also received additional interest for resistive-switching 
memoried^. As for other applications, even though the hardness of hafnia (Hf02) is not that 
high for it to be considered as a superhard materiaP, it still attracts attention as a potential 
candidate for hard oxide-based material^. Unlike carbides or nitrides, oxides are more sta¬ 
ble in the oxygen atmosphere at high temperature, which is valuable for many applications. 
Many oxide ceramics, especially those involving transition metals, are promising for applica¬ 
tion as hard coatings, since metal d electrons and strong bonds dehne their remarkable me¬ 
chanical properties (high hardness, good chemical resistance, high tensile strength, and good 
fracture toughness)!^. As compared to most transition metal oxide ceramics, hafnium oxide 
ceramics exhibit enhanced mechanical properties (higher fracture toughness) and structural 
stability (low thermal conductivity). Here, we want to explore all possible stable compounds 
of Hf-0 system at pressure up to 120 GPa. 

Under ambient temperature, experiments^ indicated that pure Hf is stable in the a- 
phase (hexagonal close-packed structure, space group: PQ^/mmc) and transforms to u- 
phase (hexagonal structure, space group: PQ/mmm) at 46-58 GPa and then to /3-Hf (body 
centered cubic, space group: ImSm) at 71.1 GPa-78.4 GPa. Our GGA calculated results 
indicate the transition pressures: a-phase —)■ cn-phase at 49 GPa and cj-phase —)■ /9-Hf 
at 70 GPa, which are in accord with above experimental results. It has been suggested 
that the solubility of oxygen in the octahedral interstitial sites of a-Hf (hcp-Hf) can be 
as high as 20 at.9c!^, while solubility of oxygen in /3-Hf (bcc-Hf) is only 3 at.9cP. Several 
experimentaP^^ and theoretical studieP^^ have investigated the interstitial oxygen in hcp- 
Hf. Now it is well established that three stoichiometric compositions HfgO, HfaO and Hf20 
can be formed with increasing occupation of the octahedral-interstitial positions in hcp-Hf by 
oxygen atoms. Hf 203 was theoretically predicted to form upon increasing the concentration 
of oxygen vacancies in monoclinic HfO^P. 

The phase sequence of Hf02 at ambient temperature with increasing pressure is: bad- 
deleyite (monoclinic, space group: P2i/c) —?• orthorhombic I (orthorhombic, space group: 
Pbca, 01) —>■ orthorhombic H (orthorhombic, space group: Pnma, OH|™. Orthorhombic 
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0II-Hf02 with experimentally reported hardness between 6-13 GPa!^^ has been speculated 
to be much harder than the low-pressure phases (baddeleyite and 0 I-Hf 02 ) because of its 
comparatively high bulk modulu^^®. 

In this study, we systematically investigate the structure and stability of Hf-0 compounds 
up to a pressure of 120 GPa by the hrst-principles evolutionary algorithm USPEX. Several 
new stoichiometries in the Hf-0 system have been predicted under high pressure. Further¬ 
more, we verify the dynamical and mechanical stability of these new high-pressure phases 
at 0 GPa by calculating their phonons and elastic constants. To better understand the 
correlations between hardness and O content, we estimate the hardness of these phases at 0 
GPa using Ghen’s hardness modeP^. Quenchable high-pressure phases often possess superior 
mechanically properties, and we indeed hnd novel hafnium oxides with unusual mechanical 
properties. 


II. COMPUTATIONAL METHODOLOGY 

Searching the stable high-pressure structures in Hf-0 system was done using hrst- 
principles evolutionary algorithm (EA) as implemented in the USPEX codd^^^ com¬ 
bined with ab initio structure relaxations using density functional theory (DFT) with 
the PBE-GGA functionaP^, as implemented in the VASP packagd^. In our work, variable- 
composition structure searcheP^ for the Hf-0 system with up to 20 atoms in the unit cell 
were performed at 0 GPa, 10 GPa, 20 GPa, 30 GPa, 40 GPa, 50 GPa, 60 GPa, 70 GPa, 
80 GPa, 90 GPa, 100 GPa, 110 GPa and 120 GPa. The initial generation of structures was 
produced randomly using space group symmetry, each subsequent generation was obtained 
by variation operators including heredity (40%), lattice mutation (20%), random (20%) 
and transmutation (20%). The electron-ion interaction was described by the projector- 
augmented wave (PAW) pseudopotentialP^, with and 2s^2p^ shells treated as 

valence for Hf and O, respectively. The generalized gradient approximation (GGA) in 
the Perdew-Burke-Ernzerhof formP^ was utilized for describing exchange-correlation effects. 
The plane-wave energy cutoff was chosen as 600 eV and P-centered uniform /c-meshes with 
resolution 2t[ x 0.06 A-i were used to sample the Brillouin zone, resulting in excellent 
convergence. Phonon dispersions were calculated using the hnite-displacement method with 
the Phonopy codd^. 
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III. RESULTS AND DISCUSSIONS 


A. Crystal structure prediction for the Hf-O system 



FIG. 1. Convex hull diagrams for the Hf-0 system at (a) 0 GPa, (b) 20 GPa, (c) 30 GPa, (d) 
40 GPa, (e) 50 GPa, and (f) 110 GPa, respectively. Solid squares denote stable phases while open 
squares represent metastable phases. 


Thermodynamic convex hull, which dehnes stable compounds, is based on the free en¬ 
ergies (at T = 0 K, enthalpies) of the compounds and pure elements in their stable forms. 
The high-pressure convex hull and pressure-composition phase diagram of the Hf-0 system 
are depicted in Fig. [T] and Fig. respectively. Besides the three well-known phases of 
Hf 02 and three suboxides (i^S-HfeO, i^Sc-HfsO and P31m-Hf20), our structure searches 
found hitherto unknown compounds with new stoichiometries, including Hf 502 , Hf 302 , HfO 


5 














1 

0.8 

o 

? 0.6 

o 

+ 

>*- 

^ 0.4 
O 

0.2 


0 

0 20 40 60 80 100 120 

Pressure (GPa) 

FIG. 2. Pressure-composition phase diagram of the Hf-0 system. 
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and HfOs- Note that a new high-pressnre phase of Hf20 (denoted as Pnnm-Hf20) was also 
fonnd by onr searches. Recent worhp! indicated that Hfi 205 is a stable componnd at low 
temperatnre, disproportionate above 220 K, therefore it is not expected to be observed ex¬ 
perimentally. Our work indicates that Hfi 205 is actually stable only in the pressure range 
from 8 GPa to 37 GPa. 


Our calculation confirms that Hf 203 proposed by Xud^can exist as a metastable phase (it 
is dynamically and mechanically stable at 0 GPa), and shows that it should be a stable phase 
in the pressure range 20-23 GPa. The predicted transition from monoclinic-Hf02 to OI-Hf02 
occurs at 10 GPa, which is coincides with experimental observation^. The transition from 
OI-Hf02 to OII-Hf02 occurs at 18 GPa, which is lower than the experimental result 30-37 
GPgpnSl good agreement with other theoretical estimates of 17 GPaP. Furthermore, 

our calculated result shows that OII-Hf02 is stable up to at least 120 GPa, which agrees with 
previous experimental worlP. According to our predictions, only baddeleyite-type Hf02 and 
RS-HfgO are stable at 0 GPa, in contrast with the Zr-0 system (ZrgO, ZrsO, Zr20, ZrO and 
Zr 02 are stable at 0 GPa;^. 
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FIG. 3. (Color online) Equation of state of Hf20. Our calculations were fit to a third-order 
Birch-Murnaghan equation of state to find Bq and Bq. 

In order to study the ordering of interstitial oxygen atoms in hcp-HfO^;, Hirabayashi et 

al.^ used electron, neutron and X-ray diffraction to analyze single crystals containing 13.4 

at % O and 15.8 at % O and found two types of interstitial superstructures: IlfOi_ and 
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HfOi , below 600 K. The space group of HfOi_ reported by Hirabayashi^ is i?3, which is 
identical to our hndings. The space group of HfOi_,_ is P31c in Hirabayashi’ experimenlP^. 
At 0 K and 0 GPa, our results produce three energetically competitive phases for HfsO and 
their ordering by energy is iJSc-HfsO (-10.075 eV/atom) < P 3 lc-Hf 30 (-10.072 eV/atom) 
< F 6322 -Hf 30 (-10.069 eV/atom). Therefore, one can note that PSlc-HfsO (P31c-Zr30 
type), exhibits very close but higher energy than i^Sc-HfsO at 0 GPa and 0 K. In order 
to consider the effects of temperature, quasi-harmonic free-energy of i? 3 c-Hf 30 and P31c- 
Hf 30 were calculated using the Phonopy cod^. The results indicate that free energy of 
F31c-Hf30 decreases faster than that of i? 3 c-Hf 30 with temperature, enabling P 3 lc-Hf 30 
to become more stable than i? 3 c-Hf 30 at 1000 K, thus explaining experimental result. 


Hf20 undergoes a trigonal-to-orthorhombic phase transition at 58 GPa. The crystal 
structure of the new high-pressure phase Pnnm-Hf20 is of anti-GaGl 2 -type. The Birch- 
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FIG. 4. (Color online) Calculated phonon dispersion curves for the (a) /mm 2 -Hf 502 at 0 and 60 
GPa (b) Pnnm-Hf20 at 0 and 80 GPa (c) G 2 /m-Hf 302 at 0 and 90 GPa (d) P62m-HfO at 0 and 
50 GPa (e) P 4 m 2 -Hf 203 at 0 and 20 GPa (f) Pnnm-HfOs at 0 and 110 GPa. The solid black and 
dashed blue lines represent the results at zero and high pressures, respectively. 


Murnaghan equation of stat^^ was used to fit the compressional behavior of the predicted 
Hf20 phases ( Fig. The third-order Birch-Murnaghan EOS is given as 


P{V} 




5 -| 

3 


1 + 4(^0 - 4 ) 



( 1 ) 


Three parameters are used to describe the EOS: the volume at 0 GPa (Vq), the bulk 
modulus at 0 GPa (Bq), and the first pressure derivative of the the bulk modulus at 0 GPa 
(Bq). Most materials have 3 < Bg < The Bg of P31m-Hf20 and Pnnm-Hf20 is 4.0 

and 3.8, respectively. 


B. Structure character in Hf-O compounds 

Table 1^ lists the detailed crystallographic data of Imm2- Hf 502 , Pnnm- Hf20, C2/m- 
Hf 302 , P62m-HfO, G 2 /m-Hf 203 and Pnnm-Hf 03 compounds at 0 GPa. The dynamical 
stabilities of all the new phases are checked by calculating phonon dispersion. As shown in 
Figj^ except for Hf 03 , no imaginary phonon frequencies are found in the whole Brillouin 
zone at both ambient and high pressure, which means that they are dynamically stable and 
probably quenchable to ambient pressure. In contrast, HfOs is stable only at high pressure, 
but at 0 GPa shows total dynamical instability and most likely decomposes. The special 
electronic structure of Hf 03 will be discussed below. The weighted average lengths of Hf-Hf 
and Hf-0 bonds in Hf-0 compounds are plotted in Fig. 
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TABLE I. Structural parameters of Imm2- Hf 502 , Pnnm- Hf20, C' 2 /m-Hf 302 , P62m-HfO, Cljm- 
Hf 203 and Pnnm-iifOs at 0 GPa. 

Compound Space group Enthalpy of formation Lattice constants Wyckoff positions x y z 


(eV/atom) (A) 


HfsOa 

Imm2 

-1.52 

a=14.455 

Hf 4c 

0.711 

0.50 

0.566 




c=3.141 

Hf 2b 

0.00 

0.50 

0.098 





Hf 4c 

0.097 

0.00 

0.594 




c=5.082 

0 4c 

0.645 

0.00 

0.818 

Hf20 

Pnnm 

-1.76 

0=5.092 

Hf 4g 

0.263 

0.341 

0.50 




6=5.723 

0 2c 

0.00 

0.50 

0.00 




c=3.175 





Hf302 

C'2/m 

-2.04 

0=11.967 

Hf 4i 

0.625 

0.50 

0.007 




6=3.131 

Hf 4i 

0.465 

0.00 

0.346 




c=11.198 

Hf 4i 

0.286 

0.00 

0.676 




d = 99.67° 

0 4i 

0.378 

0.50 

0.607 





0 4i 

0.787 

0.00 

0.192 

HfO 

P62m 

-2.60 

0=5.230 

Hf lb 

0.00 

0.00 

0.50 




c=3.187 

Hf 2c 

0.667 

0.333 

0.00 





0 3g 

0.00 

0.592 

0.50 

Hf203 

PAm2 

-3.11 

0=3.137 

Hf 2g 

0.00 

0.50 

0.744 




c=5.638 

0 2g 

0.00 

0.50 

0.135 





0 Ic 

0.50 

0.50 

0.50 

Hf03 

Pnnm 

-2.22 

0=5.554 

Hf 4c 

0.246 

0.110 

0.250 




6=6.457 

0 4c 

0.359 

0.426 

0.250 




c=3.307 

0 4c 

0.025 

0.339 

0.750 


Structurally, hafnium oxides can be divided into four groups: suboxides with oxygen in¬ 
terstitials in hcp-Hf (HfeO, Hf 30 , Hfi 205 and P31m-Hf20); other suboxides (Hf 502 , Pnnm- 
Hf20 and HfO); normal oxides (Hf 203 , Hf02) and oxide peroxide (HfOs). The octahedral 
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FIG. 5. (Color online) Average bond lengths in Hf-0 compounds at 0 GPa. 


sites of hep hafnium metal are depicted in Figj^ Oxygen atoms prefer to occupy these octa¬ 
hedral sites and form ordered structures iiS-HfgO, iiSc-HfaO, i? 3 -Hfi 205 and P31m-Hf20, 
as shown in Fig. (a) (b) (c) and (d), where Hf atom sites are omitted. The polyhedral 
representation of these structures is shown in Fig. (e) (f) (g) and (h). Anti-CaCl 2 -type 
(Fnnm) structure of Hf20 can also be represented as an hcp-sublattice (distorted) of Hf 
atoms, where half of octahedral voids are occupied by O atoms. The structure of Hf 302 can 
be considered to be defective because each layer lacks some Hf atoms to form a Hf-graphene 
layer. These vacancies are responsible for low values of the mechanical properties of Hf 302 . 


Similar with ZrO, the structure of HfO contains Hf-graphene layers stacked on top of 
each other (Zr-Zr distances within the layer are 3.01 A, and between the layers 3.18 A), 
as illustrated in Fig. 8(b), as well as additional Hf and O atoms. The structure can be 
represented as ca-phase of Hf, intercalated with oxygen atoms. This structure, therefore, 
is built by a 3D-framework of short and strong Hf-0 bonds, reinforced by rather strong 
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FIG. 6. (Color online) Octahedral voids in hep hafnium. 


Hf-Hf bonds. The former lead to high hardness, the latter may improve toughness due to 
semimetallic behavior. P 4 m 2 -Hf 203 , which was hrstly proposed by Xu^^, has 8 -fold and 
6 -fold coordination of Hf atoms, as shown in Fig[^ 


Fnnm-HfOa becomes stable at pressures above to 110 GPa. This high-pressure phase 
originally derives from oxygen atom dissolving in both octahedral and tetrahedral voids of 
a heavily distorted hep-Hf, as shown in Fig 9(a). However, due to short distances between 
tetrahedral voids in the hep structures, some O atoms form pairs and as a result HfOs 
simultaneously contains oxide and peroxide [ 0 - 0 ]^“ anions, and can be described as 
’’oxide peroxide”. The 0-0 bond length in HfOs is 1.44 Aat 110 GPa, which is a little smaller 
than the 0-0 bond length in peroxide [0-0]^“ ion with 1.47 A®' at ambient conditions. It 
seems that peroxides and oxide peroxides (e.g. AI 4 O 7 and AIO 2 ) become stabilized in many 
systems under pressurd^. 


11 







































O Hfatom 
• O atom 
o Octahedral sites 



(c)P31m-HfcO 



(d)P31m-Hf20 



(a)R^Hfi205 (b)R3-Hfi205 (e)R3-Hf60 (f)R3-Hf60 (g)R3c-Hf30 (h)R3c-Hf30 

FIG. 7. (Color online) Oxygen sublattice representation (arrangement of oxygen atoms in the 
octahedral interstitial sites) and polyhedral representation of (a)&:(b) i? 3 -Hfi 205 (c)&:(d) P31m- 
Hf20 (e)&:(f) fiS-HfgO (g)&(h) iiSc-HfaO. Oxygen-centered octahedra and oxygen vacancies are 
shown in red and pink polyhedra, respectively. Oxygen sublattice representations (b, d, f, h) show 
only oxygen atoms (filled circles) and vacancies (open circles). 



(c) C2/m-Hf302 (d) P4m2-Hf203 


FIG. 8. (Color online) Crystal structures of (a) /TOTO 2 -Hf 502 (b) P62m-HfO (c) P 4 m 2 -Hf 203 
and (d) P 4 m 2 -Hf 203 . 0-centered octahedra and 0-centered tetrahedra are shown in red and blue 
polyhedra, respectively. Large spheres-Hf atoms; small spheres-0 atoms. 
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FIG. 9. (Color online) Crystal structure of (a) Pnnm-HfOa (b) Pnnm-Hf20; (c) Pnnm-HfOs (d) 
Pnnm-Hf20. 


C. Mechanical properties of Hf-O compounds 

Previous suggested that dense high-pressure phase 0 II-Hf 02 is quenchable 

to ambient conditions and has a high bulk modulus, and might be superhard (H > 40 GPa). 
However recent study^ reported that the hardness of 0 II-Hf 02 is well below 40 GPa and 
therefore this phase is not superhard. Interestingly, our systematic results not only conhrm 
known hardness of Hf 02 polymorphs: H(OII) < H(MI) < H(OI), but also suggest that HfO 
has the highest hardness among all hafnium oxides, see Fig. [I0|(f). In addition, Pnnm-Hf 20 
and /mm 2 -Hf 502 also exhibit higher hardness than other Hf-0 compounds, as shown in Tab 
[TTj The hardness of Hf-0 compounds does not monotonically change with O content, but a 
maximum at HfO. 

The calculated modulus P, shear modulus (7, Young’s modulus E, Poisson’s ratio v, and 
hardness of all stable Hf-O compounds are depicted in Tabl^Il| and Figj^(for comparison, 
the elastic data of the high-pressure phase Pnnm-HfOs are reported at 0 GPa although it 
is unstable at 0 GPa.) From Figj^we can conclude that the high O content in the crystal 
does not guarantee high hardness of Hf-0 compounds and the structure plays an important 
role in determining mechanical properties as we discussed above. The Vickers hardness was 
calculated according to Ghen’s modeP^ 


Hv = 2*{k‘^* - 3 


( 2 ) 
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TABLE II. Calculated bulk modulus B, shear modulus G, Young’s modulus E, Poisson’s ratio v 
and hardness of Hf-0 compounds, compared with literature data for Hf 02 at 0 GPa. All properties 
are in GPa (except dimensionless G/B and v ). 


Compound 


Space group P 

Bh 

Gh 

E 

G/B 

V 


HfeO 

This work 

R3 

0 

129.2 

72.8 

183.8 

0.56 

0.26 

9.55 

Hf 30 

This work 

R3c 

0 

150.3 

78.8 

201.3 

0.52 

0.28 

9.1 

Hf502 

This work 

/mm 2 

0 

150.0 

95.3 

235.9 

0.64 

0.24 

13.9 

Hfi205 

This work 

R3 

0 

163.3 

94.5 

237.7 

0.58 

0.26 

12.1 

Hf 20 

This work 

P31m 

0 

175.2 

103.1 

258.6 

0.59 

0.25 

13.2 

Hf 20 

This work 

Pnnm 

0 

173.0 

110.3 

272.9 

0.64 

0.23 

15.5 

Hf302 

This work 

C2/m 

0 

154.2 

75.9 

195.6 

0.49 

0.29 

8.0 

HfO 

This work 

P62m 

0 

210.7 

128.1 

319.5 

0.61 

0.25 

16.1 

Hf203 

This work 

P4m2 

0 

243.9 

127.1 

324.8 

0.52 

0.28 

12.9 

HfOs 

This work 

P2ilc 

0 

203.6 

99.2 

256.1 

0.49 

0.29 

9.7 


Experiment!^ 


0 






9.9 

HfOs 

This work 

Phca 

0 

225.9 

115.8 

296.6 

0.51 

0.28 

11.7 

HfOs 

This work 

Pnma 

0 

226.3 

93.8 

247.3 

0.41 

0.31 

7.2 


ExperimentP3 


0 






6-13 

HfOg 

This work 

Pnnm 

0 

171.1 

73.6 

193.0 

0.43 

0.31 

6.2 


We calculated the elastic anisotropy of five special phases: P62m-HfO, Pnnm-Hf20, 


/mm 2 -Hf 502 , P31m-Hf20 and P 4 m 2 -Hf 203 . As shown in Fig. 11, all of these hve phases 
exhibit a moderate amount of anisotropy of Young’s modulus. The directional dependence 
of the Young’s modulus for hexagonal, orthorhombic, trigonal and tetragonal crystals can 
be calculated as: 


;— — Sii(l — ^ 3 )^ + S 33/3 + (2si3 + ■ 544 )^ 3(1 — ^ 3 ) (3) 


1 


E 


ortho 


+ 522 (^ 2 )^ + -533 (/s)^ + ^ 2 ^ 3(2523 + 544 ) + 1111(2313 + S 55 ) + Ill1(2si2 + Sge) (4) 
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O / (Hf+0) ratio 


FIG. 10. (Color online) Compositional dependence of the computed mechanical properties of Hf-0 
compounds. The blue open circle represents Pnnm-Hf20; red open triangle represents 0I-Hf02 ; 
green open square represents 0 II-Hf 02 . 


—— — (1 — / 3 )sii + / 3 S 33 + /^(l — ll){2si3 + S 44 ) + 2/2^3(3^i — ^2)'^14; (5) 

-^tri 

— - = Sii(/4 + I 2 ) + 533^3 + (2Si2 + + (2Si 3 + S44)/3(l — (6) 

^tetra 

where sn, S 12 , etc., are the elastic compliance constants and Zi, Z 2 , Z 3 are the direction cosines 
of a particnlar crystallographic orientation to coordinate axes xi, X 2 and X 3 , respectively. 
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FIG. 11. (Color online) Orientational dependence of Young’s moduli (in GPa) of (a) P62m-HfO 
(b) Pnnm-HfaO (c) /mm 2 -Hf 502 (d) PSlm-HfaO and (e) P 4 m 2 -Hf 203 . 

D. Electronic structure of Hf-O compounds 


Fig. 13 shows band structures of Hf-0 compounds at 0 GPa (including phases stable at 
both zero and high pressure). Total and partial densities of states (DOS) are presented in 


Fig. i^a-HfgO, iiSc-HfsO, Imm2- Hf 502 , i?3- Hfi 205 , P31m-Hf20, Pnnm- Hf 20 and 
G 2 /m-Hf 302 are predicted to be metallic with a sizable density of states at the Fermi level. 
The DOSs of i^S-HfgO, i? 3 c-Hf 30 , Imm2- Hf 502 , i?3- Hfi 205 , P31m-Hf20, Pnnm- HfaO 
and C2/m-Y{{^02 below Ei? are mainly due to Hf-d and 0-p orbitals, and the interactions 
between the Hf-d orbitals are responsible for metallicity. 



FIG. 12. (Color online) ELF isosurface (ELF = 0.62) for Hf 03 . Blue and red atoms represent 
oxide 0^~ and peroxide [0-0]^“ ions, respectively. 
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FIG. 13. (Color online) Band structures of hafnium oxides at 0 GPa. The Fermi energy is set to 
zero. 


Unlike regular metals, semimetals possess both electronic and hole conduction, which can 
be seen in the band structure as overlap of the partially vacant valence band top and occupied 
conduction band bottom located at different points in the Brillouin Zon^^l. Moreover, the 
upper limit of electron and hole density for a semimetal should below 10^^ cm“^. The carrier 
concentrations of the most common semimetals, for example Bi, Sb and As are 3x10^^ cm“^, 
5x10^® cm“^, and 2x10^° cm”®®, respectively. Semimetallic behavior of P62m-HfO and 
F 4 m 2 -Hf 203 can be reflected in the calculated band structures, as shown in Fig. Oh-i) 
and in the DOS diagrams (Fig. mii-i )), showing very few states at the Fermi level. For 
the band structure of P62m-HfO, there are small electron and hole pockets at F and M, 
respectively, but there are no band crossings between the lowest unoccupied bands and the 
highest occupied bands. Both band edges are mostly derived from Hf 5d states. In the case of 
F 4 m 2 -Hf 203 , partially occupied valence band top and conduction band bottom correspond 
to different high symmetry points R and Z, respectively. The electron and hole densities can 
be estimated by integrating the occupation numbers in the 3D Brillouin zone. The electron 
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FIG. 14. (Color online) The normalized (per electron) total (TDOS) and partial densities of 
states (PDOS) of hafnium oxides at 0 GPa. The Fermi energy is set to zero. 


and hole densities of HfO are both 1.1x10^^ cm by integrating their occnpation of the 


bine and green bands shown in Fig. 13 h), respectively. For P 4 m 2 -Hf 203 , the electron and 
hole densities are both 2.1x10^^ cm“^, which is very close to 1.8x10^^ cm“^ obtained in 
Xne’s workP. 

The DFT band gaps of P2i/c-Hf02, P6ca-Hf02, Pnma-Hf02 and Pnnm-HfOa are 4.01 
eV, 4.18 eV, 3.36 eV and 1.92 eV, respectively, and the highest occnpied states are all de¬ 


rived mainly from 0-p orbitals, as shown in Fig. 14 j, k, 1 and m). Therefore, according 


to their electronic character, Hf-0 componnds can be divided into three types: metallic, 
inclnding i^S-HfgO, i^Sc-HfsO, Imm2- Hf 502 , R3- Hfi 205 , P31m-Hf20, Pnnm- Hf20 and 
C' 2 /m-Hf 302 ; semimetallic, inclnding P62m-HfO and P 4 m 2 -Hf 203 ; insnlating or semicon- 
dncting, inclnding P2i/c-Hf02, F’&ca-Hf02, Pnma-IM.02 and Pnnm-]M.O^. Electron local¬ 
ization fnnction (ELF) clearly reveals special featrnre of HfOs, the coexistence of oxide 
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and peroxide [0-0]^ anions (Fig. 12). The peroxide is responsible for gap states, which 


signihcantly rednce the electronic band gap of Hf02 (Fig. 14 (m)). To obtain fnrther insight, 
we applied the Atoms in Molecnles (AIM) theory developed by Badei®. Bader charges are 
+2.5 for Hf, -0.68 for peroxide anion and -1.16 for oxide anion in HfOa at 110 GPa, which 
shows a signihcantly ionic character of bonding. 


IV. CONCLUSIONS 


We have systematically predicted stable componnds and crystal structnres in the Hf-0 
system at pressnres np to 120 GPa using ab initio evolutionary algorithm USPEX. Several 
new stable compounds, including /mm 2 -Hf 502 , G 2 /m-Hf 302 , P62m-HfO and Pnnm-HfOs 
are found for the hrst time. Pnnm-Hf 20 , which is the new high-pressure phase of Hf 20 , 
is also discovered. HfOs shows interesting structure, simultaneously containing oxide 0^~ 
and peroxide [0-0]^“ anions. Semimetallic properties of P62m-HfO and P 4 m 2 -Hf 203 are 
demonstrated through their band structures, as well as low densities of conduction electrons 
and holes. Our results demonstrate that Hf 302 is more ductile than other Hf-0 compounds, 
and the hardest compound is HfO instead of 0II-Hf02 . The superior mechanical properties 
of P62m-HfO, such as bulk modulus B, shear modulus G, Young’s modulus E and hard¬ 
ness Hy, can be attributed to the peculiar combination of strong Hf-0 and Hf-Hf bonds. 
Pnnm-Hf 20 , /mm 2 -Hf 502 , P31m-Hf20 and P 4 m 2 -Hf 203 also show excellent mechanical 
properties. Glearly, high O content is not a key factor affecting the mechanical properties 
of Hf-0 compounds. Suboxides: HfgO, Hf 30 , Hfi 205 and P31m-Hf20 based on hcp-Hf 
sublattice provide easy pathways for absorbing or desorbing oxygen. The recognition of 
the common structural features between P62m-HfO and cn-Hf gives further insight into the 
physical properties and suggests that HfO can be made as a hard semimetallic coating on 
cn-Hf substrate. Pnnm-Hf20, /mm 2 -Hf 502 , P31m-Hf20 and P 4 m 2 -Hf 203 phases in partic¬ 
ular can be quenched to ambient pressure and can be candidates for applications requiring 
mechanically strong materials. 
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